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Plasmon-X is a project based in the Frascati INFN laboratories using the Ti:Sa laser FLAME and electrons' linac 
SPARC. FLAME (Frascati Laser for Acceleration and Multidisciplinary Experiments) has a pulse with maximum 
energy Em = 7 3, maximum duration tm = 25 fs, maximum power Wm = 250 TW, wavelength A = 800 nm, and 
repetition rate v rep = 10 Hz (see, f.i.p}. The pulse duration r = 25 fs corresponds to the pulse length L z = 7.5 ^m 
(~ 10 wavelengths within the pulse). An upgrade that includes the pola rizat ion control (S, P, circular) is planned for 
the near future^. The plasma density in different Plasmon-X experiments^] ranges as n e = 0.6— 1 x 10 19 cm~ 3 , up to 
4 x 10 19 cm~ 3 . Note that n e = 10 19 cm" 3 corresponds to the plasma frequency v p ^ e = oj p ^/2tt ~ 28THz. Thus, the 
pulse duration is r = 0.7 T p , where T p is the plasma period T p = 2ir/uj p _ e = 35 fs, while the collisionless skin depth, 
d e = ~ 10.6 /im, is an order of magnitude longer than the laser wavelength (A = 0.8 fxm) and close to the 

pulse length (L z = 7.5 pxa). The laser wake field (LWF) accelerator scheme^ envisages that an electron bunch with 
the energy of 150 MeV and the transverse normalized emittance of 1 mm-mrad, whose transverse and longitudinal 
rms sizes are 5 /jm and 2.5 //m, respectively, is injected in the second bucket of the Langmuir wave excited by a Ti:Sa 
pulse delivering 7 J of energy in 30 fs. The laser pulse with initial waist size of L± = 130 pm and minimum size of 
Lj_ = 32.5 /im is guided by a matched channel profile. Plasma is 9.88 cm long and its density profile has a positive 
and varying slope with starting and ending densities of 1.5 x 10 17 cm -3 and 2.5 x 10 17 cm~ 3 . 

We derive the nonlinear equations that describe the propagation of ultrashort laser pulses in a plasma, in the 
Plasmon-X device. We consider the interaction of the high frequency electromagnetic and Langmuir waves, while the 
acoustic phenomena are disregarded. The laser electric field is so strong that the electrons achieve relativistic jitter 
velocities. The nonlinear effects come mostly from the interaction between the electromagnetic pump wave with a 
Langmuir wave, whose frequency is considerably lower than that of the electromagnetic (laser) pump. The electrons 
are regarded as cold, i.e. the phase velocity of the nonlinear modes involved are much higher than the electron 
thermal velocity, see. e.g. RefsP™H The characteristic frequency of the laser light is sufficiently high and the ions are 
essentially immobile, i.e. rij = n and Vi = 0. The component of the Ampere's law perpendicular to the direction of 
propagation of the laser beam and the Poisson's equation have the form 
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where p and j± are the charge density and the perpendicular component of the current density j, respectively, 
that are calculated as p = ^ a q a n a and j = '^2 a q a n a Vai where q a is the charge of the particle species a, and 
the hydrodynamic densities n a and velocities v a are calculated from the appropriate hydrodynamic equation. The 
electron continuity and momentum equations take the form 



dn 
dt 



+ V • (nv) = 0, 



d_ 
dt 



V p 



dA 

~dt 



+ v x Vxi 



(2) 



where, for simplicity, the subscript for electrons has been omitted and q = — e. Here p is the electron momentum, 
related with the electron velocity v through v = p/m j, is the electron rest mass, 7 = (1 + p 2 / m^c 2 )? , and c is 
the speed of light. Under the Plasmon-X conditions, the sol ution is slowly varying in the reference frame moving with 
the velocity e z u and we can use the approximation from^. Using the dimensionless quantities p — > p/m^c, v — > v/c, 
<p — > qcft/rriQC 2 , A ~ > qA/iriQC, n — > ti/uq, u — > u/c, t — > ui p , e t, f — > (u> p ^/c)(f — e z ut), where w p , e is the plasma 
frequency of stationary electrons, oj p . e = (noq 2 /nioeo)^ . The solution of the hydrodynamic equations ^ is sought in 
the almost 1-D, i.e. Vj_ <C d/dz, quasistatic regime, i.e. d/dt <C ud/dz, and u is adopted close to the speed of light, 
i.e. 1 — u <ti 1. Then, with the accuracy to the leading order, they are integrated as 



{v z - 1) n + 1 = 0, - Pz + 7 - 1 + 4> = 0, p ± + A_l = 0, A z = 0. 
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Using the definition for 7, after some straightforward algebra, we obtain the dimensionless charge- and current densities 
as n = [(<fi — l) 2 + A± + l]/2 (<fi — l) 2 and v±n = A±/(<f> — 1), which permits us to rewrite our basic equations (lib as 
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The above equations constitute a Zakharov-like description of a modulated electromagnetic wave, coupled with a 
Langmuire wave via the nonlinearities that arise from the relativistic effects. Thus, besides the standard nonrelativistic 
three- wave coupling phenomena (the Raman scattering), in the relativistic case there arises also a possibility for 
the four-wave processes, that may lead to the modulational instability, soliton formation, etc. The actual nonlinear 
dynamics of the pulse strongly depends on the physical conditions in each particular device and can not be generalized. 
We apply our Eqs. Q and |5| to the Plasmon-X conditions, for moderately focussed laser beams. With the present 

power of the laser in the Plasmon-X device, in most experimental setups we have Aj_ < 1, i.e we can take the 
pulse to have a moderate intensity and expand into series the nonlinear terms in the above equations. We seek the 

solution in the form of a modulated electromagnetic wave, viz. Aj_ = Aj_ exp{— i[uj't — k'(z + ut)]} + c.c, where 
the dimensionless frequency uj' and the dimensionless wavenumber k' are defined as uj' = uj/uj pei k' — ckjuj pe — 
d e /X, where oj, k, and A are the frequency, the wavenumber, and the wavelength of the electromagnetic laser wave, 
respectively, while uj Pyf , and d e are the electron plasma frequency and the collisionless skin depth. They satisfy the 
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whose dimensionless version has the form 



linear dispersion relation of electromagnetic waves, 

uj = \Jk 2 + 1. Here and in the following, for simplicity, we drop the primes. We adopt u to be equal to the group 
velocity of the electromagnetic wave u = dco/dk = k/u. Substituting these into the wave equation and the Poisson's 
equation, dropping the nonresonant zero- and double pump frequency terms in Eq. Q (the latter is absent for a 
circularly polarized laser wave), these are further simplified to 
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We studied the 2-D evolution of a moderately focussed pulse, typical for the Plasmon-X accelerator scheme, by the numerical 
solution of Eqs. ^ in the 2-D regime = d 2 /dx 2 . The initial condition was adopted in the form of an unchirped NLS soliton, 
modulated in the perpendicular direction by a Gausssian, Aj_ g (x, z,0) = 1\[C\ exp (— x 2 /2L 2 ) exp (i 5k z) sech(v / CT^), with 
L x — 25, Ci = 0.07, and 5k = 0.5, with (f> (x, z, 0) = 0. We followed its evolution until t = 15. During this time, the pulse 
travels approximately 4.5 cm, which is about half the length of the interaction chamber. The results are displayed in Fig. 1. 
The folding of the pancake-like pulse and the creation of a V shape was observed after t ~ 5 and the simultaneous longitudinal 
stretching due to the negative chirp. The related diminishment affected mostly the laser envelope, while the electrostatic 
potential still featured a sizable amplitude; the depth of its first minimum was more than 50% of its largest value, achieved 
immediately after the launch of the pulse. The potential minimum obtained a V-shape in the x-z plane and its perpendicular 
contraction due to collapse was rather small. 
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FIG. 1: The electrostatic potential (left) and the amplitude (right) of a moderately focussed laser pulse, typical for the 
Plasmon-X accelerator scheme. 



In conclusion, we have shown that the experimental scheme for the LWF electron acceleration, with the laser power of 
W < 250 TW, the pulse duration r > 25fs, and the spot size of 130 /jm, can be described by a nonlinear Schrodinger equation 
with a reactive nonlocal nonlinear term, that produces an oscillating electrostatic wake to the laser pulse. Such moderately 
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focussed pulses develop a chirp, that is roughly proportional to the intensity. The pulse spreads, relatively slowly, in the 
direction of propagation and its amplitude diminishes, and simultaneously a non-selfsimilar transverse collapse takes place. We 
do not expect that the latter constitutes a critical limitation for the Plasmon-X accelerator scheme. For very large amplitudes, 
~ 30 times bigger than the potential in Fig. 1, the collapse is quenched by the saturation of the nonlinear term in the wave 
equation. The presence of the chirp offers a new venue for the stabilization of the collapse even at smaller amplitudes, known 
in nonlinear optics and in Bose-Einstein condensate, and it permits one to introduce some sort of external control, such as the 
appropriate profile of the plasma density, to stabilize the structure. 



Acknowledgements. This work was partially supported by the grant 171006 of the Serbian MSE and partially by Fondo 
Affari Internazionali of INFN. 



* Electronic address: diovanov@ipb.ac.rs 

' Electronic address: renato.fedele@na.infn.it 

* Electronic address: tanjia@na.infn.it 

' Electronic address: sergio.denicola@ino.it 

1 G. Turchetti et al, Optical acceleration activity in Italy: Plasmonx, Prometheus project, 
http:/ / wwwapr.kansai.jaea.go.jp/pmrc en/org/colloquium/download/colloquiuml6-l.pdf 

2 C. Benedetti et al, NTA PtomonA, |http://www.lnf.mfn.it/rapatt/2008/06/PLASMONX.p'dr1 

3 L. A. Gizzi et al, Eur. Phys. J. Spec. Top., 175, 3 (2009). 

4 P. Sprangle, E. Esarey, and A. Ting, Phys. Rev. Lett., 64, 2011 (1990). 

5 V. I. Berezhiani and S. M. Mahajan, Phys. Rev. Lett. 73, 1837 (1994). 

6 A. Sharma, I. Kourakis, and P. K. Shukla, Phys. Rev. E 82, 016402 (2010). 



